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Abstract 

We describe Haskell implementations of interesting combi- 
natorial generation algorithms with focus on boolean func- 
tions and logic circuit representations. 

First, a complete exact combinational logic circuit syn- 
thesizer is described as a combination of catamorphisms and 
anamorphisms. 

Using pairing and impairing functions on natural number 
representations of truth tables, we derive an encoding for Bi- 
nary Decision Diagrams (BDDs) with the unique property 
that its boolean evaluation faithfully mimics its structural 
conversion to a a natural number through recursive appli- 
cation of a matching pairing function. 

We then use this result to derive ranking and unranking 
functions for BDDs and reduced BDDs. 

Finally, a generalization of the encoding techniques to 
Multi-Terminal BDDs is provided. 

The paper is organized as a self-contained literate Haskell 
program, available at http : //logic . csci .unt . edu/tare.u/ 
[research/2008/ f BDD . zip| 

Keywords exact combinational logic synthesis, binary de- 
cision diagrams, encodings of boolean functions, pair- 
ing/impairing functions, ranking/ unranking functions for 
BDDs and MTBDDs, declarative combinatorics in Haskell 

1. Introduction 

This paper is an exploration with functional programming 
tools of ranking and unranking problems on Binary Deci- 
sion Diagrams. The practical expressiveness of functional 
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programming languages (in particular Haskell) are put at test 
in the process. The paper is part of a larger effort to cover 
in a declarative programming paradigm, arguably more el- 
egantly, some fundamental combinatorial generation algo- 
rithms along the lines of (Knuth 2006). 

The paper is organized as follows: 

Sections|2]and|4]overview efficient evaluation of boolean 
formulae in Haskell using bitvectors represented as arbitrary 
length integers and Binary Decision Diagrams (BDDs). 

Section [3] describes an exact combinational circuit syn- 
thesizer. 

Section [5] discusses classic pairing and unpairing opera- 
tions and introduces new pairing/unpairing functions acting 
directly on bitlists. 

Section [6] introduces a novel BDD encoding (based on 
our unpairing functions) and discusses the surprising equiv- 
alence between boolean evaluation of BDDs and the inverse 
of our encoding, the main result of the paper. 

Section [7] describes ranking and unranking functions for 
BDDs and reduced BDDs. 

Section|8]extends our results to Multi-Terminal BDDs. 

Sections [9] and [H)| discuss related work, future work and 
conclusions. 

The code in the paper, embedded in a literate program- 
ming LaTeX file, is entirely self contained and has been 
tested under GHC 6.4.3. 

2. Evaluation of Boolean Functions with 
Bitvector Operations 

Evaluation of a boolean function can be performed one bit 
at a time as in the function if _then_else 

if _then_else _ z = z 
if _then_else 1 y _ = y 

resulting in 

> [( [x,y,z] ,if_then_else x y z) 

x^[0,l] ,y<- [0,1] ,z^[0,l]] 
[([0.0,0],0), 



([0,0,1] ,1), 
([0,1,0] ,0), 
([0,1,1] ,1), 
([1,0,0] ,0), 
([1,0,1] ,0), 
([1,1,0] ,1), 

([1,1,1] ,1)] 

Clearly, this does not take advantage of the ability of modern 
hardware to perform such operations one word a time - with 
the instant benefit of a speed-up proportional to the word 
size. An alternate representation, adapted from (Knuth 2006 ) 
uses integer encodings of 2 n bits for each boolean variable 
xq, . . . , % n -\- Bitvector operations are used to evaluate all 
value combinations at once. 

Proposition 1. Let x k be a variable for < k < n where 
n is the number of distinct variables in a boolean expression. 
Then column k of the truth table represents, as a bitstring, 
the natural number: 



Xk 



(2 2 " - l)/(2 2 



1) 



(1) 



For instance, if n = 2, the formula computes xq — 3 = 
[0,0,1,1] andxi = 5 = [0,1,0,1]. 

The following functions, working with arbitrary length 
bitstrings are used to evaluate the [0..n-l] variables Xk with 
formula [T] and map the constant 1 to the bitstring of length 
2", 111. .1: 

— the k-th, out of n bitvector boolean variables 
var_n n k = var_mn (bigone n) n k 

— the k-th, out of n boolean variables w.r.t mask 
var_mn mask n k = mask 'div' (2~ (2~ (n-k-1) )+l) 

— represents constant 1 as 11... 1 
bigone nvars = 2~2~nvars - 1 

We have used in var_n an adaptation of the efficient 
bitstring-integer encoding described in the Boolean Evalu- 
ation section of ( |Knuth|2 006 ). Intuitively, it is based on the 
idea that one can look at n variables as bitstring representa- 
tions of the n columns of the truth table. 

Variables representing such bitstring-truth tables (seen 
as projection functions) can be combined with the usual 
bitwise integer operators, to obtain new bitstring truth tables, 
encoding all possible value combinations of their arguments. 
Note that the constant is represented as while the constant 
1 is represented as 2 2 — 1, corresponding to a column in the 
truth table containing ones exclusively. 

3. Exact Combinational Circuit Synthesis 

A first application of these variable encodings is combina- 
tional circuit synthesis, known to be intractable for anything 
beyond a few input variables. Clearly, a speed-up by a factor 
proportional to the machine's wordsize matters in this case. 



3.1 Encoding the Primary Inputs 

First, let us extend the encoding to cover constants 1 and 0, 
that we will represent as "variables" n and n+1 and encode 
as vectors of n zeros or n ones (i.e. 2 2 — 1, passed as the 
precomputed parameter m to avoid costly recomputation). 



encode_var m n k 
encode_var m n k 
encode_var m n k 



k==n = m 
k==n+l = 
var_mn m n k 



Next we can precompute all the inputs knowing the num- 
ber n of primary inputs for the circuit we want to synthesize: 

init_inputs n = 

0:m:(map (encode_var m n) [0..n-l]) where 
m=bigone n 

>init_ inputs 3 

[0,15,3,5] 
>init_ inputs 3 

[0,255,15,51,85] 

Given that inputs have all distinct encodings, we can 
decode them back - this function will be needed after the 
circuit is found. 

decode_var nvars v j v=(bigone nvars) = nvars 
decode_var nvars = nvars+1 
decode_var nvars v = head 

[k | k^ [0 . .nvars-1] , (encode_var m nvars k)=v] 
where m=bigone nvars 

>map (decode_var 2) (init_inputs 2) 
[3,2,0,1] 

>map (decode_var 3) (init_inputs 3) 
[4,3,0,1,2] 

We can now connect the inputs to their future occurrences 
as leaves in the tree representing the circuit. This means 
simply finding all the functions from the set of inputs to 
the set of occurrences, represented as a list (with possibly 
repeated) values of the inputs. 

bindings us = [ [] ] 
bindings n us = 

[zs | ys^bindings (n-1) us,zs<— map (:ys) us] 

>bindings 2 [0,3,5] 

[[0,0] , [3,0] , [5,0] , [0,3] , [3,3] , 
[5,3] , [0,5] , [3,5] , [5,5]] 

For fast lookup, we place the precomputed value combina- 
tions in a list of arrays. 

generateVarMap occs vs = 

map (listArray (0,occs-l)) (bindings occs vs) 

>generateVarMap 2 [3,5] 

[array (0,1) [(0,3) , (1,3)] , 
array (0,1) [(0,5) , (1,3)] , 
array (0,1) [(0,3) , (1,5)] , 
array (0,1) [(0 ,5) , (1 , 5)] ] 



3.2 The Folds and the Unfolds 

We are ready now to generate trees with library operations 
marking internal nodes of type F and primary inputs marking 
the leaves of type V. 

data T a = V a | F a (T a) (T a) deriving (Show, Eq) 

Generating all trees is a variant of an unfold operation 
(anamorphism). 

generateT lib n = unfoldT lib n 

unf oldT 1 k = [V k] 
unfoldT lib n k = [F op 1 r 

i<-[l. .n-1] , 

1 <— unfoldT lib i k, 

r <- unfoldT lib (n-i) (k+i) , 

op<— lib] 

For later use, we will also define the dual fold operation 
(catamorphism) parameterized by a function f describing 
action on the leaves and a function g describing action on 
the internal nodes. 

f oldT _ g (V i) = g i 
foldT f g (F i 1 r) = 

f i (foldT f g 1) (foldT f g r) 

This catamorphism will be used later in the synthesis process 
for things like boolean evaluation. A simpler use would be 
to compute the size of a formula as follows: 

fsize t = foldT f g t where 
g _ = 

f _ 1 r = 1+1+r 

A first use of foldT will be to decode the constants and 
variables occurring in the result: 

decodeV nvars is i = V (decode_var nvars (is!i)) 

decodeF ixy=Fixy 

decodeResult nvars (leaf DAG, varMap,_) = 

foldT decodeF (decodeV nvars varMap) leafDAG 

The following example shows the action of the decoder: 

>decodeV 2 (array (0,1) [(0, 5) , (1 ,3) ] ) 

V 1 

>decodeV 2 (array (0,1) [ (0, 5) , (1 ,3) ] ) 1 

V 

>decodeResult 2 ((F 1 (V 0) (VI)), 

(array (0,1) [(0 , 5) , (1 ,3)] ) , 4) 
F 1 (V 1) (V 0) 

The following function uses foldT to generate a human 
readable string representation of the result (using the opname 
function given in Appendix): 

showT nvars t = foldT f g t where 

g i = 

if i<nvars 

then "x"+f(show i) 



else show (nvars+l-i) 
f i 1 r =(opname i)4+" ("4+14+" , "+fr4+") " 

> showT 2 (F 4 (V 0) (F 1 (V 1) (V 0))) 
"xor(x0,nor(xl,x0)) " 

3.3 Assembling the Circuit Synthesizer 

A Leaf-DAG generalizes an ordered tree by fusing together 
equal leaves. Leaf equality in our case means sharing a 
primary input variable or a constant. 

In the next function we build candidate Leaf-DAGs by 
combining two generators: the inputs-to-occurrences gen- 
erator generateVarMap and the expression tree generator 
generateT. Then we compute their bitstring value with a 
foldT based boolean formula evaluator. The function is pa- 
rameterized by a library of logic gates lib, the number of 
primary inputs nvars and the maximum number of leaves it 
can use maxleaves: 

buildAndEvalLeaf DAG lib nvars maxleaves = [ 
(leafDAG , varMap , 

foldT (opcode mask) (varMap!) leafDAG) 
k<— [1 . .maxleaves] , 
varMap^generateVarMap k vs , 
leafDAG ^generateT lib k 
] where 

mask=bigone nvars 
vs=init_ inputs nvars 

We are now ready to test if the candidate matches the speci- 
fication given by the truth table of n variables ttn. 

f indFirstGood lib nvars maxleaves ttn = 
head [r | r<— 

buildAndEvalLeaf DAG lib nvars maxleaves, 

testspec ttn r 
] where 

testspec spec (_,_,v) = spec=v 

> f indFirstGood [1] 2 8 1 

(F 1 (F 1 (V 0) (V 1)) (F 1 (V 2) (V 3)), 
array (0,3) [(0, 5) , (1 , 0) , (2 ,3) , (3, 0)] , 1) 

The final steps of the circuit synthesizer consist in convert- 
ing to a human readable form the successful first candidate 
(guaranteed to be minimal as they have been generated by 
increasing order of nodes). 

synthesize_f rom lib nvars maxleaves ttn = 
decodeResult nvars candidate where 
candidate^ indFirstGood lib nvars maxleaves ttn 

synthesize_with lib nvars ttn = 

synthesize_from lib nvars (bigone nvars) ttn 

— synthesizes an shows a circuit 
syn lib nvars ttn 
(show ttn)4+":"+f 

(showT nvars (synthesize_with lib nvars ttn)) 



— shows all circuits synthesized 

— for functions with nvars inputs 
synall lib nvars = 

map (syn lib nvars) [O..(bigone nvars)] 

The following example shows a minimal circuit for the 2 
variable boolean function with truth table 6 (xor) in terms 
of the library with opcodes in [0] i.e. containing only the 
operator nand. Note that codes for functions represent their 
truth tables i.e. 6 stands for [0,1,1,0]. 

> syn [0] 2 6 

"6 :nand (nand (xO.nand (xl , 1) ) , nand (xl , nand (xO, 1) ) ) " 

The following examples show circuits synthetized for 3 ar- 
gument function if-the-else in terms of a few different 
libraries. As this function is the building block of boolean 
circuit representations like Binary Decision Diagrams, hav- 
ing perfect minimal circuits for it in terms of a given library 
has clearly practical value. The reader might notice that it 
is quite unlikely to come up intuitively with some of these 
synthesized circuits. 

> syn symops 3 83 
"83:nor(nor(x2,x0) ,nor (xl ,nor (x0,0) ) ) " 

>syn asymops 3 83 

"83 : impl (impl (x2 , xO) , less (xl , impl (xO , 0) ) ) " 
>syn mixops 3 83 

"83:nand(impl(x2,x0) ,nand(xl ,x0) ) " 

> syn [3,4] 3 83 

"83 : xor (xl , less (xor (x2 , xl ) , xO) ) " 

We refer to the Appendix for a few details, related to the 
bitvector operations on various boolean functions used in the 
libraries, as well as a few tests. 

4. Binary Decision Diagrams 

We have seen that Natural Numbers in [0..2 2 — 1] can be 
used as representations of truth tables defining n-variable 
boolean functions. A binary decision diagram (BDD) ( |Bryant| 
1986) is an ordered binary tree obtained from a boolean 
function, by assigning its variables, one at a time, to (left 
branch) and 1 (right branch). 



The construction is known as Shannon expansion ( Shan- 
|non|19 93 ), and is expressed as a decomposition of a function 
in two cof actors, f[x<—0] and f[x «— 1] 



f(x) = (xAf[x^0])V(xAf[x^l}) 



(2) 



where f[x <— a] is computed by uniformly substituting a 
for x in /. Note that by using the more familiar boolean 
if-the-else function, the Shannon expansion can also be 
expressed as: 



f(x) — if x then f[x <— 0] else f[x <— 1] 



(3) 



Alternatively, we observe that the Shannon expansion can 
be directly derived from a 2" size truth table, using bitstring 
operations on encodings of its n variables. Assuming that the 



first column of a truth table corresponds to variable x, x = 
and x = 1 mask out, respectively, the upper and lower half 
of the truth table. 

Seen as an operation on bitvectors, the Shannon expan- 
sion (for a fixed number of variables) defines a bijection as- 
sociating a pair of natural numbers (the cof actors' s truth 
tables) to a natural number (the function's truth table), i.e. 
it works as a pairing function. 

5. Pairing Functions 

Definition 1 . A pairing function is a bijection f : Nat x 
Nat — > Nat. An unpairing function is a bijection g : 

Nat -> Nat x Nat 

5.1 Classic Pairing Functions 

Following Julia Robinson's notation (Robinson 1950), given 
a pairing function J, its left and right inverses K and L are 
such that 



J(K(z),L(z)) = z 



K(J(x,y))=x 



(4) 



(5) 



L(J(x,y))=y (6) 

We refer to ( Cegielski and Richard 20011) for a typical use 
in the foundations of mathematics and to ( Rosenberg 2002 ) 
for an extensive study of various pairing functions and their 
computational properties. 

Starting from Cantor's pairing function 



f(x,y) = (x + y)* (x + y+l)/2- 
and the Pepis-Kalmar-Robinson function 

f(x,y) = 2* *(2*y+l)-l 



(7) 



(8) 



bijections from Nat x Nat to Nat have been used for 
various proofs and constructions of mathematical objects 



(Pepis 


19381 |Kalmar| 


1939} 


Cegielski and Richard 


2001) 



5.2 Pairing/Unpairing operations acting directly on 
bitlists 

We will introduce here a pairing function, expressed as sim- 
ple bitlist transformations. This unusually simple pairing 
function (that we have found out recently as being the same 
as the one in defined in Steven Pigeon's PhD thesis on 
Data Compression (Pigeon 200 l|l, page 114), provides com- 



pact representations for various constructs involving ordered 
pairs. 

The function bitmerge_pair implements a bijection 
from Nat x Nat to Nat that works by splitting a num- 
ber's big endian bitstring representation into odd and even 
bits, while its inverse bitmerge_unpair blends the odd and 



even bits back together. The helper functions nat2set and 
set2nat, given in the Appendix, convert from/to natural 
numbers to sets of nonzero bit positions. 

bitmerge_pair (i,j) = 

set2nat ((evens i) +f (odds j)) where 
evens x = map (2*) (nat2set x) 
odds y = map succ (evens y) 

bitmerge_unpair n = (f xs,f ys) where 
(xs.ys) = partition even (nat2set n) 
f = set2nat . (map ('div' 2)) 

The transformation of the bitlists is shown in the follow- 
ing example with bitstrings aligned: 

>bitmerge_unpair 2008 
(60,26) 

— 2008: [0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1] 

— 60: [ 0, 1, 1, 1, 1] 

— 26: [ 0, 1, 0, 1, 1 ] 

PROPOSITION 2. The following function equivalences hold: 
bitmerge jpair o bitmerge _unpair = id (9) 

bitmerge -impair o bitmerge jpair = id (10) 

6. Pairing Functions and Encodings of 
Binary Decision Diagrams 

We will build a BDD by applying bitmerge_unpair re- 
cursively to a Natural Number tt, seen as an n-variable 2™ 
bit truth table. This results in a complete binary tree of depth 
n. As we will show later, this binary tree represents a BDD 
that returns tt when evaluated applying its boolean opera- 
tions. 

We represent a BDD in Haskell as a binary tree BT 
with constants and 1 as leaves, marked with the function 
symbol C. Internal nodes representing if -then-else de- 
cision points, marked with D, are controlled by variables, 
ordered identically in each branch, as first arguments of D. 
The two other arguments are subtrees representing the THEN 
and ELSE branches. Note that, in practice, reduced, canoni- 
cal DAG representations are used instead of binary tree rep- 
resentations. 

data BT a = C a | D a (BT a) (BT a) 
deriving (Eq, Show) 

The constructor BDD wraps together the number of vari- 
ables of a binary decision diagram and the binary tree repre- 
sentation it. 

data BDD a = BDD a (BT a) deriving (Eq, Show) 

The following functions apply bitmerge_unpair recur- 
sively, on a Natural Number tt, seen as an n-variable 2™ bit 
truth table, to build a complete binary tree of depth n, that 
we will represent using the BDD data type. 



— n=number of variables, tt=a truth table 
plain_bdd n tt = BDD n bt where 

bt=if tt<max then shf bitmerge_unpair n tt 
else error 

("plain_bdd: last arg "+f (show tt)4+ 
" should be < " +f (show max)) 
where max = 2~2~n 

— recurses to depth n, splitting tt into pairs 
shf f n tt n<l = C tt 

shf f n tt = D k (shf f k ttl) (shf f k tt2) where 
k=pred n 
(ttl,tt2)=e tt 

The following examples show the results returned by plain_bdd 
for all 2 2 truth tables associated to n variables for n = 2, 
with help from printing function print_plain given in Ap- 
pendix. 

>print_plain 2 
BDD 2 (D 1 (D (C 0) (C 0)) (DO (C 0) (C 0))) 
BDD 2 (D 1 (D (C 1) (C 0)) (DO (C 0) (C 0))) 
BDD 2 (D 1 (D (C 0) (C 0)) (DO (C 1) (C 0))) 

BDD 2 (D 1 (D (C 0) (C 1)) (D (C 1) (C 1))) 
BDD 2 (D 1 (D (C 1) (C 1)) (D (C 1) (C 1))) 

6.1 Reducing the BDDs 

The function bdcLreduce reduces a BDD by collapsing 
identical left and right subtrees, and the function bdd asso- 
ciates this reduced form to n £ Nat. 

bdd_reduce (BDD n bt) = (BDD n (reduce bt)) where 
reduce (C b) = C b 

reduce (D _ 1 r) | 1 = r = reduce 1 

reduce (D v 1 r) = D v (reduce 1) (reduce r) 

bdd n = bdd_reduce . plain_bdd n 

Note that we omit here the reduction step consisting in 
sharing common subtrees, as it is obtained easily by replac- 
ing trees with DAGs. The process is facilitated by the fact 
that our unique encoding provides a perfect hashing key for 
each subtree. 

The following examples show the results returned by bdd 
for n=2, with help from printing function print_reduced 
given in Appendix. 

>print_r educed 2 
BDD 2 (C 0) 

BDD 2 (D 1 (D (C 1) (C 0)) (C 0)) 
BDD 2 (D 1 (C 0) (D (C 1) (C 0))) 
BDD 2 (D (C 1) (C 0)) 

BDD 2 (D 1 (D (C 0) (C 1)) (C 1)) 
BDD 2 (C 1) 



6.2 From BDDs to Natural Numbers 

One can "evaluate back" the binary tree representing the 
BDD, by using the pairing function bitmerge_pair. The 
inverse of plainJbdd is implemented as follows: 



plain_inverse_bdd (BDD 
rshf bitmerge_pair bt 



bt) = 



rshf rf (C tt) = tt 

rshf rf (D _ 1 r) = rf ((rshf rf 1) , (rshf rf r) ) 



>plain_bdd 3 42 
BDD 3 
(D 2 
(D 



(D 1 



(D (C 0) 

(D (C 0) 

(D (C 1) 

(D (C 1) 



(C 0)) 

(C 0))) 

(C 1)) 

(C 0)))) 



plain_inverse_bdd it 
42 

Note however that plain_inverse_bdd does not act as an 
inverse of bdd, given that the structure of the BDD tree is 
changed by reduction. 

6.3 Boolean Evaluation of BDDs 

This rises the obvious question: how can we recover the orig- 
inal truth table from a reduced BDD? The obvious answer is: 
by evaluating it as a boolean function! The function ev de- 
scribes the BDD evaluator: 

ev (BDD n bt) = eval_with_mask (bigone n) n bt 

eval_with_mask m _ (C c) = eval_constant m c 
eval_with_mask m n (D x 1 r) = 
ite_ (var_mn m n x) 

(eval_with_mask m n 1) 

(eval_with_mask m n r) 

eval_constant = 
eval_constant m 1 = m 

The function ite_ used in eval_with_mask implements 
the boolean function if x then t else e using arbitrary 
length bitvector operations: 

ite_ x t e = ( (t 'xor' e).&.x) 'xor' e 

We will use ite as the basic building block for implement- 
ing a boolean evaluator for BDDs. 

6.4 The Equivalence 

A surprising result is that boolean evaluation and structural 
transformation with repeated application of pairing produce 
the same result, i.e. the function ev also acts as an inverse of 
bdd and plainJbdd. 

As the following example shows, boolean evaluation ev 
faithfully emulates plain_inverse_bdd, on both plain and 
reduced BDDs. 



>plain_bdd 3 42 
BDD 3 
(D 2 
(D 



(D 1 

ev it 
42 

bdd 3 42 
BDD 3 
(D 2 
(C 0) 



(D 
(D 
(D 
(D 



(C 0) 
(C 0) 
(C 1) 
(C 1) 



(C 0)) 

(C 0))) 

(C 1)) 

(C 0)))) 



(D 1 
(C 1) 

(D (C 1) 



(C 0)))) 



ev it 
42 



The main result of this subsection can now be summa- 
rized as follows: 

PROPOSITION 3. The complete binary tree of depth n, ob- 
tained by recursive applications of bitmerge_unpair on 
a truth table tt computes an (unreduced) BDD, that, when 
evaluated, returns the truth table, i.e.: 



plain .inverse Jbdd (plainJbdd ntt)) = id 



(11) 



ev n (plainMd n tt)) = id (12) 
Moreover, ev also acts as a left inverse of bdd, i.e. 



ev n (bdd ntt)) = id 



(13) 



Proof sketch: The function plainJbdd builds a binary tree 
by splitting the bitstring tt 6 [0..2™ — 1] up to depth n. 
Observe that this corresponds to the Shannon expansion 
(Shannon 1993) of the formula associated to the truth table, 
using variable order [re — 1, 0]. Observe that the effect of 
bitstring_unpair is the same as 

• the effect of var_mn m n (n-1) acting as a mask select- 
ing the left branch, and 

• the effect of its complement, acting as a mask selecting 
the right branch. 

Given that 2™ is the double of 2™~ 1 , the same invariant holds 
at each step, as the bitstring length of the truth table reduces 
to half. On the other hand, it is clear that ev reverses the 
action of both plainJbdd and bdd, as BDDs and reduced 
BDDs represent the same boolean function (Bryant 1986). 

This result can be seen as a yet another intriguing isomor- 
phism between boolean, arithmetic and symbolic computa- 
tions. 

7. Ranking and Unranking of BDDs 

One more step is needed to extend the mapping between 
BDDs with n variables to a bijective mapping from/to Nat: 



we will have to "shift towards infinity" the starting point of 
each new block of BDDs in Nat as BDDs of larger and 
larger sizes are enumerated. 

First, we need to know by how much - so we will count 
the number of boolean functions with up to n variables. 

bsum = 

bsum n | n>0 = bsuml (n-1) 
bsuml 0=2 

bsuml n | n>0 = bsuml (n-l)+ 2~2~n 

The stream of all such sums can now be generated as usuaQ 

bsums = map bsum [0 . . ] 

>genericTake 7 bsums 

[0 , 2 , 6 , 22 , 278 , 65814 , 42950331 10] 

What we are really interested into, is decomposing n into 
the distance n-m to the last bsum m smaller than n, and the 
index that generates the sum, k. 

to_bsum n = (k,n-m) where 

k=pred (head [x | [0 . . ] ,bsum x>n] ) 
m=bsum k 

Unranking of an arbitrary BDD is now easy - the index k 
determines the number of variables and n-m determines the 
rank. Together they select the right BDD with plain_bdd 
and bdd. 

nat2plain_bdd n = plain_bdd k n_m 
where (k,n_m)=to_bsum n 

nat2bdd n = bdd k n_m 
where (k,n_m)=to_bsum n 

Ranking of a BDD is even easier: we shift its rank within the 
set of BDDs with nv variables, by the value (bsum nv) that 
counts the ranks previously assigned. 

plain_bdd2nat bdd@(BDD nv _) = 

(bsum nv)+(plain_inverse_bdd bdd) 

bdd2nat bdd@(BDD nv _) = (bsum nv)+(ev bdd) 

As the following example shows, nat2plain_bdd and 
plain_bdd2nat implement inverse functions. 

>nat2plain_bdd 42 
BDD 3 
(D 2 
(D 1 
(D 
(D 
(D 1 



(C 0) 
(C 1) 



(C 1)) 
(C 0))) 
(D (C 0) (C 0)) 
(D (C 0) (C 0)))) 
>plain_bdd2nat it 
42 

1 bsums is sequence A060803 in The On-Line Encyclopedia of Integer 
Sequences, http : //www . research. att . com/~nj as/sequences 



The same applies to nat2bdd and its inverse bdd2nat. 

>nat2bdd 42 
BDD 3 
(D 2 
(D 1 

(D (C 0) (C 1)) 
(D (C 1) (C 0))) 
(C 0)) 
>bdd2nat it 
42 

We can now generate infinite streams of BDDs as follows: 
plain_bdds = map nat2plain_bdd [0. .] 

bdds = map nat2bdd [0. .] 

>genericTake 4 plain_bdds 
[ 

BDD (CO), 
BDD (C 1) , 

BDD 1 (D (C 0) (C 0)) , 
BDD 1 (D (C 1) (C 0)) 

] 

genericTake 6 bdds 
[ 

BDD (CO), 
BDD (CI), 
BDD 1 (CO), 

BDD 1 (D (C 1) (C 0)) , 
BDD 1 (D (C 0) (C 1)) , 
BDD 1 (C 1) 

] 

8. Multi-Terminal Binary Decision Diagrams 
(MTBDD) 

MTBDDs ( |Fujita et al.][l997j |Ciesinski et~aT1|2008] l are a 
natural generalization of BDDs allowing non-binary values 
as leaves. Such values are typically bitstrings representing 
the outputs of a multi-terminal boolean function, encoded as 
unsigned integers. 

We shall now describe an encoding of MTBDDs that 
can be extended to ranking/unranking functions, in a way 
similar to BDDs as shown in section [7] 

Our MTBDD data type is a binary tree like the one used for 
BDDs, parameterized by two integers m and n, indicating 
that an MTBDD represents a function from [0..n — 1] to 
[0..m — 1], or equivalently, an n-input/m-output boolean 
function. 

data MTBDD a = MTBDD a a (BT a) deriving (Show.Eq) 

The function tojntbdd creates, from a natural number 
tt representing a truth table, an MTBDD representing func- 
tions of type N -> M with M = [0..2 m - 1],N = 
[0..2 n — 1]. Similarly to a BDD, it is represented as binary 
tree of n levels, except that its leaves are in [0..2 m — 1]. 



to_mtbdd m n tt = MTBDD m n r where 
mlimit=2~m 
nlimit=2~n 

ttlimit=mlimit~nlimit 
r=if tt<ttlimit 

then (to_mtbdd_ mlimit n tt) 
else error 

("bt: last arg "4+ (show tt)+f 

" should be < " -H- (show ttlimit)) 

Given that correctness of the range of tt has been checked, 
the function to_mtbdd_ applies bitmergejunpair recur- 
sively up to depth n, where leaves in range [0. .mlimit — 1] 
are created. 

to_mtbdd_ mlimit n tt | (n<l)&&(tt<jnlimit) = C tt 
to_mtbdd_ mlimit n tt = (D k 1 r) where 

(x,y)=bitmerge_unpair tt 

k=pred n 

l=to_mtbdd_ mlimit k x 
r=to_mtbdd_ mlimit k y 

Converting back from MTBDDs to natural numbers is 
basically the same thing as for BDDs, except that assertions 
about the range of leaf data are enforced. 

from_mtbdd (MTBDD m n b) = from_mtbdd_ (2~m) n b 

from_mtbdd_ mlimit n (C tt) | (n<l)&&(tt<mlimit)=tt 
from_mtbdd_ mlimit n (D _ 1 r) = tt where 
k=pred n 

x=^f rom_mtbdd_ mlimit k 1 
y=^f rom_mtbdd_ mlimit k r 
tt=bitmerge_pair (x,y) 

The following examples show that tojntbdd and f ronuntbdd 
are indeed inverses values in [0..2™ — 1] x [0..2™ — 1]. 

>to_mtbdd 3 3 2008 
MTBDD 3 3 
(D 2 
(D 1 



(D 





(C 


2) 


(C 


D) 


(D 





(C 


2) 


(C 


1))) 


(D 1 












(D 





(C 


2) 


(C 


0)) 


(D 





(C 


1) 


(C 


1)))) 


mtbdd 


it 











2008 



>mprint (to. 


.mtbdd 2 


2) 


[0. 


3] 












MTBDD 2 2 






















(D 1 (D 


(C 


0) 


(C 


0)) 


(D 





(C 


0) 


(C 


0))) 


MTBDD 2 2 






















(D 1 (D 


(C 


1) 


(C 


0)) 


(D 





(C 


0) 


(C 


0))) 


MTBDD 2 2 






















(D 1 (D 


(C 


0) 


(C 


0)) 


(D 





(C 


1) 


(C 


0))) 


MTBDD 2 2 






















(D 1 (D 


(C 


1) 


(C 


0)) 


(D 





(C 


1) 


(C 


0))) 



9. Related work 

Pairing functions have been used for work on decision prob- 
lems as early as ( |Pepis| 1 938] |Kalmar 1939; Robinson 1950| >. 

BDDs are the dominant boolean function representa- 
tion in the field of circuit design automation ( |Meinel and| 
|Theobald| 19991 |Drechsler et al.|2004| i. 

Besides their uses in circuit design automation, MTBDDs 
have been used in model-checking and verification of arith- 
metic circuits (Fujita et a l.|1997| Ciesinski et al.|2 008). 

BDDs have also been used in a Genetic Programming 
context dSakanashi efaT1[T996l |Rothlauf et aLl|2006l [Chen] 
|et al.|2004| as a representation of evolving individuals sub- 
ject to crossovers and mutations expressed as structural 
transformations. 



10. Conclusion and Future Work 

Our new pairing/unpairing functions and their surprising 
connection to BDDs, have been the indirect result of imple- 
mentation work on a number of practical applications. Our 
initial interest has been triggered by applications of the en- 



codings to combinational circuit synthesis (Tarau and Luder- 
man|20 08). We have found them also interesting as uniform 
blocks for Genetic Programming applications. In a Genetic 
Programming context ( |Koza|199"2"t|Poli et al.[ ), the bijections 
between bitvectors/natural numbers on one side, and trees/- 
graphs representing BDDs on the other side, suggest explor- 
ing the mapping and its action on various transformations as 
a phenotype-genotype connection. Given the connection be- 
tween BDDs to boolean and finite domain constraint solvers 
it would be interesting to explore in that context, efficient 
succinct data representations derived from our BDD encod- 
ings. 
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Appendix 

To make the code in the paper fully self contained, we list 
here some auxiliary functions. 

Bitvector Boolean Operation Definitions 

type Nat=Integer 

nand_ : : Nat-^Nat^Nat^Nat 
nor_ : : Nat^Nat^Nat^Nat 
impl_ : : Nat^Nat^Nat^Nat 
less_ : : Nat^Nat^Nat^Nat 

nand_ mask x y = mask .&. (complement (x .&. y) ) 
nor_ mask x y = mask .&. (complement (x . | . y)) 
impl_ mask x y = (mask .&. (complement x)) . | . y 
less_ x y = x .&. (complement y) 

Boolean Operation Encodings and Names 

— operation codes 
opcode m = nand_ m 
opcode m 1 = nor_ m 
opcode m 2 — impl_ m 
opcode m 3 = less_ m 
opcode 4 = xor 

opcode n = error ("unexpected opcode : "4-f( show n) ) 

— operation names 
opname = "nand" 
opname 1 = "nor" 
opname 2 = "impl" 
opname 3 = "less" 
opname 4 = "xor" 

opname n = error ("no such opcode : "4-K show n) ) 
A Few Interesting Libraries 

mixops = [0,2] 
symops = [0,1] 
asymops = [2,3] 

Tests for the Circuit Synthesizer 

t0 = f indFirstGood symops 3 8 71 

tl — syn asymops 3 71 

t2 = mapM_ print (synall mixops 2) 

t3 = syn asymops 3 83 — ite 

t4 = syn symops 3 83 

t5 = syn [0. .4] 3 83 — ite with all ops 

— x xor y xor z — cpu intensive 
t6 = syn asymops 3 105 

Bit crunching functions 

This function splits a natural number in a set of natural 
numbers indicating the positions of its 1 bits in its right to 
left binary representation. 

nat2set n = nat2exps n where 
nat2exps _ = [] 
nat2exps n x = 

if (even n) then xs else (x:xs) where 
xs=nat2exps (div n 2) (succ x) 



http 



This function aggregates a set of natural numbers indicating 
positions of 1 bits into the corresponding natural number. 

set2nat ns = sum (map (2~) ns) 
I/O functions 

These functions print out the BDDs of all the 2 2 truth tables 
associated to k variables. 

print_plain k = mapM_ 

(print . (plain_bdd k) ) [O..(bigone k)] 
print_reduced k = mapM_ 

(print . (bdd k) ) [0 . . (bigone k)] 

This function applies f to a list of objects and prints the 
results on successive lines. 

mprint f = (mapM_ print) . (map f) 



